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SUMMARY 


In. this paper, the effect of nonideal square-law detection on static calibration 
for a class of Dicke radiometers is examined. It is shown that fourth-order curva- 
ture in the detection characteristic adds a nonlinear term to the linear calibration 
relationship normally ascribed to noise-injection, balanced Dicke radiometers. The 
minimum error, based on an optimum straight-line fit to the calibration ourve, is 
derived in terms of the power series coefficients describing the input-output charac- 
teristic of the detector. These coefficients can be determined by simple measure- 
ments, and detection nonlinearity is, therefore, quantitatively related to radio- 
metric measurement error. 


INTRODUCTION 

A microwave radiometer is a receiver specifically designed to measure microwave- 
radiated power P and thereby determine the radiometric brightness temperature T 
defined by b' 


rp = 

b kB (1) 


where k is Boltzmann's constant, and B is the bandwidth. Radiometers have been 
implemented in a variety of configurations (ref. 1) and are still in an evolutionary 
phase in which both theoretical and hardware innovations are being incorporated. 

This work concerns an error-producing mechanism in pulsed-noise, balanced Dicke 
radiometers and is a theoretical examination of the measurement error resulting from 
nonideal square- law detection in such a system. 

The square-law detector (S.L.D.) is present in all noise-injection radiometers 
and can be realized in many forms utilizing a variety of nonlinear elements. In this 
analysis, the S.L.D. was modeled as a memoryless, nonlinear "black box" having an 
input-output relationship given by 


V 

o 


k ; v i„ * k i v L * k ; v L * v 4 ln 


( 2 ) 


where V Q and V in are the detector output and input random-noise voltages, respec- 
tively, and represents the coefficients in the power series. This investigation 

was limited to the important case where the S.L.D. is very nearly square- law because 
of the values of the coefficients in equation (2) and a restriction on |v. I. An 
earlier investigation of this subject (ref. 2) was part of a broader discussion of a 
particular radiometer and was highly condensed and less general than this analysis. 


SYMBOLS 


B 

dBR 

E{ } 
k 



bandwidth 


» 10 log,, -i 
10 P 


r 

expected value operator 
Boltzmann's constant 


coefficient of second-order term-in power series relating v a and p 
coefficient of fourth-order term in power series relating v 0 and P, 
coefficient of first-order term in power series relating V<J and V*, 
coefficient of second-order term in power series relating v o and V. 
coefficient of third-order term in power series relating V J and vj 
coefficient of fourth-order term in power series relating v o and V, 


loss factor 


number of decades on logarithmic scale 
power 

antenna (input) power 

error associated with straight-line approximation to p ( 0 ) 

€1 

approximation to p 

d 

optimum linear approximation to p ( 0 ) 
detector input power, radiometer mode 
injecfed-noise source power 
reference power at detector input 
system noise power at detector input 
detector input power, sinusoidal excitation 
detector input resistance 
DicRe-cycle period 
antenna -(input) temperature 
temperature error corresponding to . Ap 

A 


2 



z 

Az 

o 



injected-noise source temperature 
reference termination temperature 
Byatem equivalent input noise temperature 
time 

antenna (input) random-noise voltage 
detector (input) random-noise voltage 
injected random-noise voltage 
detector output random-noise voltage 
reference-noise voltage 
system noise voltage 

detector output during first interval of Dicke cycle 
detector output during Second interval of Dicke cycle 
detector output during third interval of Dicke cycle 
value of V QX designated as "reference level" 
value of detector output voltage, sinusoidal input 
magnitude of sinusoidal test signal 

logarithmic transformation of V QX , referenced to v or 
deviation from straight line 
ratio of k 4 to k 2 
injected-noise duty cycle 

A 

value of 0 for which P a = 0 

normalized departure from linearity, reference-noise excitation 
normalized departure from linearity, sinusoidal excitation 
power-balance error on the ith Dicke cycle 
standard dev- ation of antenna noise voltage 
standard deviation of injected-noise voltage 
standard deviation of system noise voltage 


3 


ORIGINAL PAGE IS 
OF POOR QUALITY 

max maximum 

min minimum 

S.L.D, square- law detector 

Mathematical notation: 

= very nearly equal to 


Abbreviations : 

F.S. full scale 


THEORY 


Null-Balance Static-Error Derivation 


In this section, an expression is derived to relate the input (antenna) power 
and the pulsed-noise duty cycle in a pulsed noise-injection, balanced Dicke radiom- 
eter, which is shown in simplified form in figure 1. The system is a closed- loop, 
nulling servomechanism, which adds noise to the antenna noise such that their sum 
power, as measured over the antenna half-Dicke cycle, equals the power measured over 
the reference half-Dicke cycle. This analysis deals specifically with the case of a 
single noise pulse per Dicke cycle? however, the results are applicable to a pulse- 
rate modulation format. Furthermore, the target is assumed fixed, such that the 
input noise process can be considered as stationary and ergodic. 


Detector waveforms typifying the radiometer model of figure 1 
state are shown in figure 2. The system is designed to force 
average value of a sequence of random variables, to approach zero, 
able £ ^ is defined as follows: 


in an equilibrium 
the finite-time 
The random vari- 


t 4 2 f 

L T M i-1 


(i-1 )T+T/2 


)T 


V (t) dt 
o 


-C, 


V_(t) dt 


>T+T/2 


(3) 


< 


This is accomplished by the error-nulling nature of the servomechanism, which drives 
the duty cycle of the injected noise & to be that value necessary to make 
approach zero. Equation (3) can be rewritten so that each integral is restricted to 
a time interval in which the detector input consists of only one combination of the 
four relevant noise processes: input noise V (t), injected noise V (t>, system 

noise V g (t), and reference noise V y (t). Thus, n 
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at 


♦§r 

T Mi-i 


(i-1/2)T 
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where the -subscripts on V Q (t) distinguish the three distinct intervals of a Dicke 
cycle, as depicted in figure- (2) » To the extent that the servo system drives Ej 
to zero and, as previously stated, the input noise processes, are all stationary such 
that time and ensemble averages can be interchanged, the "power -balance condition" 
requires that 


Applying this constraint to- e quation ( 4) and integrating yields 


P E{V Qi (t)} + (1 - j3) E{V o2 (t)} - E{V o3 (t)} - 0 (6 ) 


During the first time interval, the detector input consists of the sum of the 
antenna noise, injected noise, and system noise, and the first term of equation (6) 
may be written as 



(7) 


following equation (2). The random processes V a , v n , and V g are assumed to be 
Gaussian, zero-mean, and statistically independent. Further reduction of equa- 
tion (7) is based on the application of 

e(v^| =( * 3 * 5 * " 1 (for j even)! 

1 1 (for j odd)j <8) 

from reference 3, and the fact that E{pq} = E{p} • E{q> if p and q are indepen- 
dent random variables. Expanding equation (7) and applying these results leads to 


(S =(V o1 I 




2 

2 J n 


+ kid 2 
2 s 


+ 6k*o 2 a 2 + 6k*a 2 o 2 + 
4 a n 4 a s 


+ 3k Jo + 
4 a 

.. . 2 2-, 
uk'o a J 
4 ns' 


3k'a^ 

4 n 


3k;<r 4 
4 s 


(9) 


This can be expressed as 




S ( k 2 P » * k 2 P n * k 2 P s 

+ 6k P P + 6k P P 
4 a n 4 a s 


* 3k„pf <■ 


+ 6k P P ) 
4 n a ' 


2 

+ 3k A P c 

4 s 
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I. 

V-' 


I' 



since k' ■ k^R and k^ - k^/R 2 , where R is the input resistance of the detector* 
li'cation of similar arguments to the ramaininn two- eXDected—va nnnraH nno 


The app 
yields 


E{V o2 } " k 2 P a + k ? P «, + 3k a p ! + 3k ^ p2 + P 

oj z a 2 s 4a - 4s 4as 


( 11 ) 


E{V o 3 } * k 2 P r + k 2 P s + 3k 4 P r + 3k 4^. t 6k 4 P r P s 


and the power-balance equation (eq, (4)) can be written as 


( 12 ) 


p ( k 2?a + k 2 P n + k 2 P s + 3k 4 P a + 3k 4 P n + 3k 4 P s + 6k 4 P a P n + 6k 4 P a P s 
+ 6k 4 P n P s) + (1 - P ) ( k 2 P a + k 2 P s + 3k 4 P a + 3k 4 P I + 6k 4 p a P s ) 


“ k 2 P r “ k 2 P s - 3 Vr ' 3k 4 P s ~ 6k 4 P r P s = 0 


(13) 


Upon normalizing by k 2 and defining a = k 4 /k 2 , equation (13) reduces to 

2. 2 
3aP + (1 + 6agP + 6aP )P + pp + 3a8P + 6a8P P 
3 n sa n r n ns 


- P - 3aP - 6aP P =0 
t r r s 


( 14 ) 


This quadratic equation in P a defines the relationship between P a and (3 i 

resulting from the "power-balance 1 * condition* In a perfect square-law detector, a ^ 

equals zero and equation (14) simplifies to a linear relationship between P 
and-0 a 


P a * p r - P p n (15) 


The major question to be answered by this paper is* What error in measuring P a 
results when a linear relationship between P a and p is assumed, as in equa- 
tion (15), when in fact, a is actually nonzero, albeit very small, and equa- * 

tion (14) more accurately applies to the situation? This is not a trivial ques- 
tion, as the usual assumption in the radiometry field is that the detector is l 

ideal, and therefore, a linear calibration relationship is justified (ref. 1, * 

p. 402). - t 

I • 

i 

A representation for P a (0) is sought which explicitly expresses any nonlinear- | 

ity between P a and f) when a is nonzero. The most direct approach to obtaining 
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such an expression for P a (3) is to-solve equation (14) for P a with the under- 
standing. that the solution will be a function of fl. Application of the quadratio 
formula to equation ( 1 4 > yields (see appendix A) 


p a “ 6aV (1 + 6a{3p n + _ 6ap s > + +_6otP r + 6aP s ) 


1 + 


36a 2 3(0 ~ OP 2 
n 


1 / 2 ' 


(16) 


(1 + 6aP + 6aP ) 
r s 


which is an exact solution only if 3 equals zero or unity, in which cases the radi- 
cal- reduces to unity,, and equation (16) yields- exactly the same value for P a as 
equation (1$). The range of 3 Is dependent on the value of P . For- a- = 0, the 
range of 3 is n 


0 < 3 < P /P 
r n 


(for P > p > 0) 
r a 


07) 


A simplified form of equation (16) that is good for all values within the range 
of 3 can be obtained by approximation. Since a has been restricted to very small 
values, it is justifiable to approximate equation (16) as 


p a “ gj<^<1 + 6 <*& p n + 6aP g ) + (1 + 6ap f + 6aP g ) 


1 + 


13a 2 3(3 - 1)P 2 
n 


( 18 ) 


(1 + 6aP -k 6ap )‘ 
r s 


which simplifies to 


3aP 2 \ 3aP 2 

2 L 2 

1-+~6<tt> + 6aP / 1 +- 6ap + 6aP 

r s/ r s 


Equation (191 shows that the fourth-order term in equation (2) has resulted in 
a slope change relative to the ideal situation where a « o given by equation (15), 
This -can be - accommodated in the calibration process), however, nonzero values of a 
also produced a nonlinear 3 2 term which-cannot be removed by a straight-line 
calibration. 

Some feeling for the physical significance of these results can be obtained by 
examining the behavior of equation (19) for three ranges of P n » P n < p r , p n * p r , 
and P n > P t , with o either zero or nonzero. The first case, P n < P r , has limited 
practical importance in radiometr-y, but it is included here for completeness. 
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Figure 3 shews P a versus p for the above-mentioned conditions. The following 
points should be noted* 

1. Hie linear (ct * 0) and nonlinear (a > 0 and a < 0) curves- intersect at 

P =» 0 -and p » 1. (See eq. (IS).) 

2. Hie practical range- of p is- bounded by zero and the lesser of approximately 

p r/ p n or 1 ' since- p fl oannot be negative, and p cannot exceed unity. 

3. In the practical operating region where P n > S r , increasingly positive 

values of a make p Q , the value of p for which P a = 0, increasingly 
less than- Pj-/P n ,. and increasingly negative values- of a make P Q 
increasingly .greater than P r /P n . (See eq. (C7).) 

4. Increasing a increases the spread between the linear and nonlinear cases. 

Now that an expression is available which -explicitly expresses P a as a func- 
tion of p, it is possible to determine the error arising: from a straight-line cali- 
bration. Hie "best" straight-line fit will be defined as that which results in the 
smallest peak error between itself and equation (19) over the region from p = 0 
to P 88 P 0 ' which corresponds to a p a range of to zero. Ihis optimum linear 

fit to -P a ( p) is designated as P a ^, ®nd the quadratic error function 


AP = P ( B ) - P , 
a a H al - 


( 20 ) 


is to be minimized in the sense stated above. It is shown in appendix B that this 
optimum linear fit to P a (P> results in a peak value for AP a that is one-eighth 
of the value attained by the p 2 term in equation (19) when p attains its maxi- 
mum value of p c , at which point equals zero. It is necessary to determine p Q 

in order to evaluate AP a , Hiis is accon^.lished in appendix C by solving for the 
root of equation (19), which is shown in equation (C7) to be 


0 


o 



3a(p - P ) 
n r 


1 + 3ap + 6aP 

n rj 


The minimum value of Ap a as described above and in appendix B can now be evaluated 
using equations (T9), (B14), and (C7), and is 

3aP 2 3a(P n - P r ) -| 2 

iP a,peaklmin T 8(V + 6op + 6aP ) 1 1 + 3aP + 6aP * 21 

r s L n sj 

which essentially equals 3ap|/8 when the degree of nonlinearity is slight. This 
is the principal result of this section and establishes a loWef bound on the radio- 
metric measurement error caused by nonzero Values of a. It should be emphasized 
that the actual error based on the usual ,K two-point“ calibration technique will 
almost certainly exceed this minimum value. 
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Detector- Nonlinearity Characterization 

Evaluation of aquation C2V) requires- a knowledge of P c and cq P r is a known 
or easily measured parameter. 'Phis section .discusses two methods for determining a. 
The first method 1 b an "in-circuit" method, suitable for use with an existing radiom- 
eter; the second method is based on sinusoidal excitation. 

In the in-circuit method, the Dicka switch is "locked" in the reference mode, 
and a calibrated step attenuator of loss factor L is placed at the input of the 
detector. The input-output characteristic is then determined by using the random 
signals V r + V g as excitation. The average value of the detector output voltage 
is given by equation (12) as 


V o3 “ e { V o 3 U) } " Y* (P r + V + 3 V 2(P r + V 2 


( 22 ) 


In figure 4, V 0 j is plotted as a function, ©f the input power L(P r + P s ). The 
normalized departure from linearity is defined as the ratio of the deviation of v 0 3 
from the linear term k 2 L(P r + P s ) | L=1 to the linear term, and oan be expressed as 


AV o3 

k 2 (P r * V 


3 V p r * V 

k 2 (P r * V 


3a(P r + P s ) 


(23) 


from which a is found to be 


3(P + P ) 

r s 


(24) 


Thus, a can be defined using the measured quantity 6-, and the known or easily 
measured total input power in the reference mode. 

Determination of a using a sinusoidal input signal follows the same general 
procedure. The mean output voltage of the detector, described by equation (2), 
when driven by =* V x sin wt (w is the radian frequency), can be found to be 


ox 


k’V 2 3k !V 4 
2 x 4 x 

2 8 


which can be written in terms of power,, since kl ^ k 2 /R, kl £ k./R 2 , and 
P x ■ V X /2R. Thus, i 


(25) 



3k, P 


k.P + 
2 x 


4 x 


(26) 
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from which the norma lized-dcpartura from linearity (see fig, — 54 can ba written- as 


3k .P 3aP 
, 4 x x 

6 2 “ 2k P “ 2 

2 x 


and a found to be 


Of. st- 


26, 

i 

3P - 


(27i 


(28! 


i 



It is new possible to numerically evaluate AP a> , as given by equation (21), 
for some typical radiometer parameter values. Suppose, for example, that 
T r = T g a 308 K, that equation (28) is evaluated at a sinusoidal power level of 
? x = P 8 + P r . an ^ i‘ lat ^2 is 0.01, which approaches the limits of analog observ- 
ability. Then, from the approximation to equation (21), 



(29 


The AT a associated with this AP a is, from equations (1) and (29), 


AT = 0.001 25T - 0.385 K 

a min r 


(30 


and it is clear that the detector must be extremely true square-law if radiometric 
errors of less than 1 K are demanded and a linear relationship between P a and p 
is assumed. Alternately, since the effect is predictable, it could be included by 
using a nonlinear calibration relationship such as equation Cl 91 » This approach 
becomes more attractive if the possibility exists that other nonideal aspects of the 
system hardware could produce nonlinearity between and p. Admitting the pres- 

ence of a- p 2 term in P a (P) could significantly improve calibration accuracy. 

As a final consideration, the consequences of following the usual laboratory 
procedure of logarithmically plotting the data used to define a (ref. 4) will be— 
examined. The resulting compression of the observed variable-signif ieantly masks 
nonlinearity. To show this, equation (26) is written such that the log reference is 

the value of v when p e p , or 
ox x r 


V k-P + (3/2)k.P 2 
ox 2-x ' 4 x 

5 „r ♦ u/2>k,p? 


(P x /P r )[l + (3<xP x /2)] 

1 + ( 3aP /2) 
r 


(31 


Taking 10 


A 

z « 10 log 


l°3l0 


ox 


10 V 


or 
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ef both aides of equation (31) yields 


*> 10 log + 

r 


10 log 




( 3ap \ 
10 10g ioV +- 2 *) 


or 


z =» dBR - 10 


log 


10 




+ 10 log 




( 32 ) 


(33) 


which is sketched in- figure 6. Now, the coefficient ratio a is evaluated in terms 
of the transformed variable z. Prom figure 6, the normalized deviation from a 
straight line at full scale (F.S.) is 


Az A 10 ^lot 1 + (3a V 2)] 

F.S. “ 10N 


P =P 
x 


r 


log 


10 




(34) 


when N is the number of decades on the vertical scale. Writing a two-term 
Maclaurin’s series approximation for the right-hand side of equation (34) in 
terms of aP r yields 

Az . 3 log 10 6 3aP r 

■— - — — as ctP = ■■ 

F.S. • 2N r r 2N In 10 


(35) 


For comparison purposes, consider the same data plotted on a linear scale. When 
P x = P r , the normalized deviation frc i a straight line was shown to be 


6 


2 


■JS 



(36), 


Comparison of equations (35) and (36) reveals that loqarithmic plotting has reduced 
the observable deviation from a straight line by the factor N In 10. In other 
words, the radiometer error associated with a given Az/F.S. has increased by the 
factor N In 10. For the example cited earlier, based on a 62 of 0,01, the AT a 
of 0.385 K would increase to 2.66 Kl Thus, logarithmic plotting should be avoided 
when attempting to define a. 
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CONCLUSIONS 


The results of the examination of the effect of nonideal square-law detection on 
the static calibration of noise-injection, balanced Dicke radiometers may be summa- 
rized as followat 

1. Fourth-order curvature in the square-law detector was shown to produce a 

nonlinear relationship between the pulsed duty cycle $ and the measured 
antenna temperature, 

2. the minimum peak-measurement error with respect to an optimum straight-line 

fit was related to a ratio of coefficients in a power series describing the 
input-output characteristic of the square-law detector. 

3. The coefficient ratio was shown to be obtainable from simple laboratory 

measurements . 

4^. The minimum calibration nonlinearity resulting from square- law detection 

which deviates from ideal by only 1 percent was s hown to b e on the order of 
0.5 K. 


5. The accuracy penalty incurred by the use of logarithmically plotted data to 
find the coefficient ratio was assessed and shown to be significant. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
November 23, 1983 
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DERIVATION OP NONLINEAR CALIBRATION EQUATION 

This. .appendix includes the steps followed in the derivation of equation (19), 
for which the quadratic solution roots are, by inspection, 

P 6aj)P + 6ctP ) ± (1 + 36a V p*; + 3&«V + 12a6P n + 12aP g 

a - -6a u n 3 

+- 72a 2 6P P - 12aBP - 36a 2 6P 2 - 72a 2 0P P + 12aP 
ns n 9 

2 „ 2 ^ o 0/ 2 l < 

+ 36a P + 72a P P J J 
r us 

After simplification and regrouping, equation (Al) may be written as 

P - ^ — [-( 1 + 6a6P + 6aP ) + [l + l2a(P + P ) + 36a 2 (P r + P g ) 2 
a 6a ns to 

+ 36a 2 6 ( 6 - 1)P 2 ] 1/2 } 


and simplified to 



+ 6a&P n + 6aP g ) ± 


(1 + 6ap + 6aP ) 1 
r s 


+ 


36a 2 6(6 - 1)P 2 
n 

2 

(1 + 6aP + 6aP ) 
r s 



Because of the stipulated smallness of a (corresponding to only slight devia^ 
tion from true square-law in the detector), the term in brackets is approximated by 
the first two terms of a Maclaurin’s series to yield 


^a “ 6a V (1 + 6aBP n + 6oP s ) * + 6aP r + 6aP s* 


1 + 


18a 2 6(6 - DP 2 


(A 


(1 + 6aP + 6aP ) 
IT S 


Only the positive square root i 3 significant and produces realistic values for P a . 
Thus, equation (A3) simplifies to 


P - P, - 6P 
a r n 


1 + 3aP + 6aP + 6aP 
n r s 


3ap‘ 


1 + 6aP + 6aP 

IT 9 


1 + 6aP + 6aP 
r s 


(? 


which is equation (19) in the text in a slightly different form. 
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An approximation differing only slightly from equation (A4) can be obtained by 
writing a Maclaurin's series for equation (14) with the coefficients obtained by 
implicitly differentiating equation (14) and evaluating the derivatives at 0 = 0, 
Shis procedure yields an approximation for P a that differs from equation (A4) only 
in the coefficient of the 0^ term# which-is 


d 2 P 
1 a 

21 d3 2 


0=0 


3aP*(1 + 3aP + 6aP + 6aP ){1 - 3oP e 6aP + 6aP ) 
n n r s n r s__ 

1 + 6aP +■ 6aP 
r s 


(A5) 


and which can be expressed as 


,2 

1 d P a 
21 d0 2 


0=0 


3aP 


1 + 


1 + 6aP + 6aP 


3 OP 


1 * 6aP + 6ap 
r s 


(A6) 


to emphasize the extreme closeness of the two approximations. 
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determination op best straight-line pit to a quadratic function 

ON A FINITE INTERVAL 1 

For a quadratic equation 
2 

y = ax — + bx + c 

the-conditions are derived that cause the linear approximation te y, 

Y = Ax + B 
L 

to be the best fit to equation (B1 ) in the sense that 

4y U* 

is minimized over the interval from zero to x m . The quadratic equation 
Ay = y - * ax 2 + {b - A)x + (c - B) 

can have a maximum of two real roots for which Ay « 0 . Thus, over the interval 
2 ero to x m , the roots are positioned so that 

Ay(0) = Ay(x ) ■ -Ay{x ) 
m o 

as shown in the sketch below* 


(B1 ) 

(B2) 

(B3) 

(B4) 

from 

(B5) 
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Applying the first set of equalities in equation (B5) to equation (B4) yields 


c 



+ (b 


A )x + c 
m 


B 


which reduces to 


(B6) 


A = b + ax m (B7) 

and defines the slope of Y T . 

Jj 

In order to find the intercept of Y^# it is necessary to first find x Q , the 
flex point over. -the interval from zero to x ra . Prom equation (B4), the derivative 




2ax + (b - A) 


(B8) 


must- equal zero at x 0 . Thus, 

A - b 

x o " ~ir~ <«) 

Combining this with the result stated in equation (b 7) yields 


m 


x = 
o 2 


(BIO) 


Thus, the value of Ay at x Q is 


Ay(x ) * a -p + (b - b - ax )~ + (c - B) 
04 mz 


(B11) 


which by specification of equation (B5) must equal 
B yields 


-y(0), or -(c - B) . solving for 


ax 


B = c - 


m 


(B1 2 ) 
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by evaluating equatio^ M) ^ei ^ by eq J att J n < B2 > can- be found 

' x o' or x m . Evaluating at zero yields 


Ay min “ l Ay(0) I - o - B 


Combining this with equation (B12) gives the desired result 


(B13) 
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DETERMINATION OP 0 Q 

In- this appendix, the value of 0 for which P a , equation (19), is zero is 
determined. Equation (19) is first put in standard quadratic form, 


. <■ 3oP ♦ 6«P t 4-6«P a 1 ♦ 6oP * 6®, 

* [ 5ST )K * —2 fr * “ 


and the quadratic formula applied tq_yield 


(l + 3ctP + 6aP + 


1 + 3aP + 6ap + 6ap 


1 2a(1 + 6oP + 6aP ) ' ' 

r s_ 


Removing common factors and expanding gives 


1 + 3oP + 6aP + 6aP 

0 1 * + 9a 2 P 2 * 36a 2 p 2 ♦ 36a 2 p 2 

o 6aP 6aP ^ nr s 

n n 

2 2 2 

+ 6aP + 12aP + 1 2oiP + 36a p P +• 36a P p + 72a P P - 1 2aP 
n r s nr ns rs r 


- 72a 2 p 2 - 72a 2 p P ) 1>/2 
r r s' 


which can be simplified and written as 


1 + 3ap + 6aP + 6aP 
n r 


1 * eirto + 3aP n + 6aP s ) 2 + 36a 2 P r (P n - P r >]‘ 


This can be put in a form suitable for approximating by removing the large factor 
from the term in brackets; thus, 


1 + 3ap + 6aP + 6aP 1 + 3ap + 6aP 


36a 2 p (P - P ) “l 1 ^ 2 


(1 + 3ap + 6aP ) 
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Figure 2.-* Input and output waveforms at detector for full Dicke 
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Figure 4.- Graphical determination of nonlinearity 
coefficient using reference-noise excitation. 



Figure 5.- Graphical determination of nonlinearity 
coefficient using sinusoidal excitation. 






